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Abstract. 

We present the results of recent observations of phase-dependent variations 
in brightness designed to characterize the atmospheres of hot Jupiters. In par- 
ticular, we focus on recent observations of the transiting planet HD 189733b 
at 8 /xm using the Spitzer Space Telescope, which allow us to determine the 
efficiency of the day-night circulation on this planet and estimate the longi- 
tudinal positions of hot and cold regions in the atmosphere. We discuss the 
implications of these observations in the context of two other successful detec- 
tions of more sparsely-sampled phase variations for the non-transiting systems 
V And b and HD 179949b, which imply a potential diversity in the properties of 
the atmospheres of hot Jupiters. Lastly, we highlight several upcoming Spitzer 
observations that will extend this sample to additional wavelengths and more 
transiting systems in the near future. 



1. Introduction 



There are currently more than 20 known transiting planetary systems, of which 
the majority are gas-giant planets orbiting extremely close (<0.05 A.U.) to their 
parent star^. These planets, known as "hot Jupiters", receive >10,000 times 
more radiation from their stars than Jupiter does froin the s un, heating them 
to temperatures as high as 2000 K ( Harrington et all l2007l ). The picture is 
further complicated by the fact that most of these planets are expected to be 
tidally locked, with permanent day and night sides. As a result, the equilibrium 
compositions and circulation patterns in these atmospheres are expected to differ 
significantly from those of the gas-giant planets in the solar system. 

It is not surprising, t hen, that siniple n i odels for the atmospheres of these 
planets (see, for exainple, ISeager et al.ll2005l : iBarman et al.l [20051 : iFortney et al 



planets isee. tor exainpie. I5e 
2006bl : iBurrows et al.ll2007bl ^ 



can sometimes differ significantly in both their as- 
sumptions and the conclusions that they reach; this is particularly true of the cir- 
culation models discussed in ^2]] Fortunately, it is possible to test the predictions 
of these models using observations of transiting systems. From the wavelength- 
dependent depth of the transit, when the planet passes in front of the star, 
we can search for features in the transmission spectrum of the planet near the 
day-night t erminator (Charb onneau et al- 2002: Vid al-Madiar et al. 2003 , 200J; 
Knutson et al.l l2007al : IBarmaneTaLl l2007i : lTinetti et al.1 l2007l : lEhrenreich et al.l 
20071 1. Similarly, by measuring the depth of the secondary eclipse, when the 
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sion from the day side of the planet ('Charbonneau et al.' '2005': 'Deming et al. 


2005 


20061. 2007: Demorv et al. 


2007: Grillmair ct al. 


2007: Richardson et al. 


2007 


Harrington et al. 20071: Knutson et al. 


2007blJcl). Lastly, by measurine: the 



changes in brightness over a significant fraction of a planet's orbit we can di- 
rectly constr ain the longitudirial temperature distribu t ion across the planet 's 
atmosphere (|Harrington et al.ll2006l : ICowan et alll2007l : iKnutson et al.n2007bl ') . 
This last type of observation is particularly informative for hot Jupiters, which 
may have extreme changes in temperature between the day and night sides of 
the planet. It is also a difficult observation to make in practice as it requires a 
very high precision over time scales of several days; this is why all of the suc- 
cessful detections to date have been made from space. Although the shape of 
the ingress and egress during secondary eclipse also cont ain information about 
the brightness distri bution on the day side of the planet ( Williams et al.l [20061 : 
Rauscher et al. 20071 ) , these variations are smaller and take place on shorter time 



scales, and there have not been any clear detections of this effect. 

In this paper we focus on observations of phase variations, discussing the 
predictions of the current generation of atmospheric circulation models, the re- 
sults of the three successful measurements of phase variations to date, and up- 
coming plans for more such observations in the near future. 



2. Predictions from Atmospheric Circulation Models 



As mentioned above, the extreme environments and slow rotation rates of the 
(presumably) tidally-locked hot Jupiters mean that their atmospheric dynamics 
are exp ected to differ significantly from those of the planets in our own solar 



Burkert et all l2005l: 



svstem (IShowman fc Gmllotl l2002l: ICho et al.1 [2003. [20061. 
Cooper Showmanll2005l . l2006l : lLangton Laughlinll2007l : 
20071 ). One of the most fundamental questions is what fraction, if any, of the 
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energy absorbed by the perpetually-illuminated day side is transferred around 
to the night side of the planet. The answer to this question depends on the rela- 
tive sizes of the radiative and advective time scales, quantities which are poorly 
constrained at the moment and may vary from planet to planet. The relative 
temperature difference that we observe between the day and night sides also 
depends on the wavelength of the observations and the opacity of the planet's 
atmosphere at that wavelength. For wavelengths where the opacity is low the 
effective photosphere of the planet is located deep in the atmosphere, where the 
pre ssures and temperatu res are correspondingly higher. In their dynamical mod- 
els Cooper &i Showman p005i ) show that circulation is more efficient at these 
higher pressures, as the advective time scale decreases relative to the radiative 
time scale (see Fig. [1]). Thus we would expect that observations of the same 
planet at different wavelengths would show varying temperature differences be- 
tween the day and night sides, depending on how deep into the atmosphere we 
are looking at each wavelength. 

The second observationally testable prediction made by circulation models 
has to do with the location of hotter and cooler regions on the planet as a 
function of longitude. In the most basic picture the hottest region is located 
at the substellar point, which has the highest levels of incident radiation, and 
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Figure 1. Atmospheric circulation model for HD 209458b from 
ICooper &: Showman! (|2005D showing the temperature distribution (grayscale) 
and winds (arrows) at three different pressures, centered on the substellar 
point. Panel (a) is at 2.5 mbar, (b) is at 220 mbar, and (c) is at 19.6 bars. 
This model predicts that circulation will be more efficient at depth, as demon- 
strated by the decreasing temperature differential and increasingly shifted 
position of the hot spot on the dayside in the lower panels. 
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Figure 2. Left panel : 24 /zm phase variation for v Andromedae b from 
[Harrington et al.l (|2006f ). with best fit phase curves overplotted (solid line is 
with no phase shift, dotte d line has a phase s hift). Right panel: 8 phase 
curve for HD 179949 from Cow an et all (|2007) . The two bottom panels show 
the 8 fim fluxes for HD 179949 and a comparison star used to correct for 
time-dependent instrument effects. The top panel shows the corrected light 
curve at for HD 179949 as a function of orbital phase with best-fit sinusoid 
overplotted. For both systems, the large amplitude of the observed varia- 
tion relative to the predicted planet-star flux ratio for some reasonable set 
of assumptions about the planet's radius, orbital inclination, and effective 
temperature is consistent with relatively inefficient circulation and a large 
day-night temperature difference. 



the coolest spot is located at the antistellar point. In atmospheres where the 
advective time scale is comparable to the radiative time scale the positions of 
these hot and cool regions are shifted in the direction of the prevailing winds. 
By determine the orbital phases where the planet reaches its maximum and 
minimum brightness, it is possible to constrain the longitudinal positions of 
these hotter and colder regions. 



3. Observing Phase Variations: Two Methods 

To date there are two observing strategies that have been pursued in the quest to 
measure the phase variations for various planets. In the first strategy, the same 
system is observed for short periods of time (usually ~ 30 minutes) at intervals 
ranging from several hours to several days. In the second method the system is 
observed continuously over half an orbit, on the order of ~ 30 hours or so. The 
advantage of this second method is that it is easier to correct for instrumental 
effects and the resulting data has a much higher cadence and correspondingly 
higher precision; however it is also much more resource-intensive. 

3.1. Method 1: A Series of Widely Spaced, Short Visits 

Harrington et all ( 20061 ) made the first detection of phase variations for an ex- 



trasolar planet by using this method to observe v Andromedae b at 24 fim using 
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the Spitzer Space Telescope (see left panel of Fig. Because this is not a tran- 
siting system we do not know the radius of the planet (although we can make an 
educated guess based on the mass), the inclination of the pla net's orbit, or the 
total f lux from the day side of the planet. As a result, although [Harrington et al" " 



( 20061 ) reported a relatively large variation in flux with orbital phase, the actual 
temperature difference between the day and night sides is not as well-constrained 
as it would be for a transiting system. If one makes reasonable assumptions, 
however, the size of the observed variation implies a large day-night temperature 
difference and correspondingly inefficient circulation between the day and night 
sides^ 

Cowan et al.l ( 20071 ) made similar observations of three other systems at 



8 iim using Spitzer, including HD 209458, HD 179949, and 51 Peg, of which 
only one (HD 209458) is a transiting system. They report a detection for one of 
the non-transiting systems, HD 179949, (see right panel of Fig. [2|) and conclude 
that the size of the observed variation implies a large day-night temperature 
difference. They also calculate upper limits on the size of the phase variations 
in the other two systems, but these upper limits place relatively weak constraints 
on the amount of circulation in the atmospheres of the two planets. 



3.2. Method 2: Continuous Observations Over Half an Orbit 



iKnutson et al. ( 2007bl ) observe the transiting planet system HD 189733 con- 



tinuously over half an orbit (33 hours) at 8 fiui using Spitzer. We began our 
observations before the start of the transit, when the night side of the planet is 
visible, and end after the secondary eclipse when the day side of the planet has 
rotated into view. Because we observe both the transit and secondary eclipse 
we are able to derive new values for the radius of the planet (from the depth 
of the transit) and the total flux from the day side (from the depth of the 
secondary eclipse), and we know we are viewing the system in an effectively 
edge-on configuration. This allows us to translate the relative increase in flux 
over the observed part of the orbit into a longitudinal temperature distribution 
across the surface of the planet, as described below. We n ote that the pres- 
ence of spots on the star, which is known to be quite active ( Winn et al. 20071: 



Henry Sz Wind 120071) . and correcti ons for detector effects (ICharbonneau et al 



2005 : Morales-Calderon et al.l [20061 ) both contribute additional uncertainties in 



th e interpretation of this measurement; both effects are discussed in more detail 



m 



Knutson et aD (|2007bl 1. 



We measure a total increase in flux of 0.12 it 0.02% over the period of our 
observations (see Fig. [3]) . If we scale this value by the depth of the secondary 
eclipse, this means that the minimum hemisphere-integrated flux is 64 it 7% of 
the maximum hemisphere-integrated flux. Converting this to brightness temper- 
atures, we find that the minimum hemisphere-averaged brightness temperature 
is 973 lb 33 K, and the maximum hemisphere-averaged brightness temperature is 
1212 lb 11 K. The uncertainties in the minimum brightness temperature, which 
approximately corresponds to the night side temperature, are larger because 
this part of the planet was visible at the beginning of the observations when the 
uncertainties in the correction for the detector ramp are largest. 

We also find that the minimum brightness occurs 6.7 zb 0.4 hours after the 
center of the transit, while the maximum brightness occurs 2.3 ib 0.8 hours be- 
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Figure 3. L e ft pan el: Observed phase variation for HD 189733b from 
iKnutson et al.l (|2007b[ ). with transit and secondary echpse visible. The stel- 
lar flux as measured at the center of the secondary eclipse is normalized to 
unity (dashed line), and the data is binned every 500 points (200 s). a and b 
show the same data, but in b the y axis is expanded to show the scale of the 
variation. The solid line is the phase curve for the best-fit model shown on 
the right. Right panel: Brightness estimates for twelve longitudinal strips on 
the surface of HD 189733b. The brightness estimates are given as the ratio of 
the flux from an individual slice viewed face-on to the total flux of the star. 



fore the center of the secondary eclipse. These shifts imply that the hottest and 
coolest points are not located at the substellar and antistellar points, but are 
instead shifted to the east and to the west, respectively. We can refine our esti- 
mates of the locations of these spots by fitting the observed phase variation with 
a model made of twelve longitudinal strips of varying brightness (Fig. [3]); the re- 
sulting brightness distribution implies that the hot spot is shifted approximately 
30 degrees east of the substellar point, while the cool spot is shifted 30 degrees 
west of the antistellar point. The shif ted hot spot is consistent with the west-to- 
east winds predicted by rnost models (iShowman &: Guillod 2002 Cho et al. 20031. 



20061 : iBurkert et al.1 [20051: ICooper k Showmanll2005l . l20od : iLangton &: Laughlin 



200?! : iDobbs-Dixon LinI l2007l 1. but these models also predict that the cool 



spot should also be shifted to the east. 



Diversity of Planetary Atmospheres There are several conclusions we can draw 
from these data. First, unlike v And b and HD 179949b, we find that for 
HD 189733 the circulation between the day and night sides of the planet must 
be relatively efficient in order to produce such a small temperature difference 
between the day and night sides. Although it is possible that some of the 
difference between v And b and HD 189733b may be due to the fact that the 
V And b observations were at a longer wavelength (24 fj,m instead of 8 jum), it 
is unlikely that this difference alone would produce such a large change in the 
amplitude of the observed phase variation. Furthermore, this cannot account 
for HD 179949b's large phase variation, as this system was observed at the same 
wavelength as HD 189733b. 
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In light of the recent detectio n of a atmosp h eric te mperature inversion 
on the dayside of HD 209458b by iKnutson et al.l ( 2007d ). it is more plausi- 
ble to think that the atmosphere of HD 189733b has fundamentally different 
properties than those of v And b and HD 179949b. There are some initial 
indications that the presence of a temperature inversion may be correlated 
with the level of irra.diation a nd /or the surface gr avity of the planet in ques- 
tion (jBurrows et alj l2007al lbl: iFortnev et alj l2007l ). as evidenced by the fact 
that TrES-1 and HP 189 7 33b do not appear t o have temperature inv ersions 
( Charbonneau et al.1 l2005l : iDeming et alJ 12006: iRichardson et al. 2007), while 



HD 209458b and HP 149026b probablv do (iGrillmair et alJl2007l: IKnutson et al. 
2007bl : iBurrows et "alll2007al : iFortnev et alll2006al : [Harrington et alJl2007l ). r; And b 



and HD 179949b both have higher levels of incident flux than HD 209458b; thus 
it would not be unreasonable to think that both of these planets may have tem- 
perature inversions. Such inversions would naturally lead to a lar ger phase vari 



ation, as the incident flux is absorbed higher in the atmosphere (jFortnev et al, 
20071 1. 



4. Conclusions and Future Directions 

The above dis cussion highlights the significant advances that have been made in 
the year since Harrington et all (j2006) reported the first detection of the phase 



variation of an extrasolar planet. However, the current sample is limited, and 
there is clearly much to be gained from observations of additional transiting 
systems, for which the planetary radius, orbital inclination, and dayside flux 
are all well-constrained. There are a number of exciting programs scheduled 
for the upcoming Spitzer observing cycle which will address this issue. First, 
we have a program whi ch will observe HD 189 733b continuously at 24 //m over 
the same half orbit as IKnutson et al.l ( 2007bl ) . which should probe a different 



region of the atmosphere than the previous 8 fim observations. As a part of 
this same program we will also observe HD 209458b over half an orbit at both 
8 and 24 /xm, which will allow us to directly compare the wavelength-dependent 
properties of the atmospheric circulation for these two planets. 

To date all of the observations we have described have examined planets 
with effectively circular orbits; however, there is another promising set of up- 
coming observations by G. Laughlin and G. Bakos that will observe two highly 
eccentric planets (HD 80606b and HAT-P-2b) continuously over ~ 30 hours cen- 
tered on periastron passage. Although only one of these planets, HAT-P-2b, is 
transiting, both observations should detect the rapid increase in brightness as 
the planetary atmosphere is flash-heated during this passage, providing a direct 
measurement of the radiative time constant. 

Of the two methods described above, it is the more time-intensive, contin- 
uous observations of transiting systems that have provided the most detailed 
information about circulation in the atmospheres of hot Jupiters to date. How- 
ever, the detections of phase variations for v And b and HD 179949b highlight the 
advantages of more sparsely-sampled observations in facilitating a much wider- 
ranging survey of a large number of planetary systems. These two approaches 
are inherently complementary, and we expect that both will make valuable con- 
tributions to our understanding of these unusual planets in the near future. 
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